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Abstract 

We establish rigourously the scaling properties of the Lanczos process applied to an 
arbitrary extensive Many-Body System which is carried to convergence n — > oo and the 
thermodynamic limit N — > oo taken. In this limit the solution for the limiting Lanczos 
coefficients are found exactly and generally through two equivalent sets of equations, 
given initial knowledge of the exact cumulant generating function. The measure and the 
Orthogonal Polynomial System associated with the Lanczos process in this regime are also 
given explicitly. Some important representations of these Lanczos functions are given, 
including Taylor series expansions, and theorems controlling their general properties are 
proven. 



PACS: 05.30.-d, 11.15.Tk, 71.10.Pm, 75.10.Jm 



E-mail: nsw@physics .unimelb . edu.au 



1 



I. Introduction 



The Lanczos Algorithm is one of the few reliable and general methods for computing the 
ground state and excited state properties of strongly interacting quantum Many-body Sys- 
tems. It has been traditionally employed as a numerical technique on small finite systems, 
with attendant round-off error problems, although the main obstacle to its further develop- 
ment has been the rapid growth of the number of basis states with system size. The reader 
is referred to a review of the applications of this method |Q] in strongly correlated electron 
problems. In this work we examine the Lanczos process in the context of the extensive quan- 
tum Many-Body Systems, where it is employed entirely in an exact manner and where the 
thermodynamic limit is taken. So in complete contrast to the traditional use of the Lanczos 
algorithm - we completely circumvent the issues of loss of orthogonality due to round-off 
errors and the inability to approach the thermodynamic limit because of the requirement to 
construct a full basis on the cluster. The systems we have in mind are those with an infinite 
number of degrees of freedom, yet are extensive, in that all total averages of any physical 
quantity scale linearly with the numbers of degrees of freedom however quantified. These 
would include all condensed matter systems with sufficiently local interactions (the precise 
conditions need to be clarified, but it is clear which specific systems obey extensivity) and 
Quantum Field Theories, with the proviso that the spectrum is bounded below (in some 
cases there is also an upper bound too). 

After noting some of the advantageous features of the algorithm in general we discuss the 
scaling behaviour of the Lanczos Process as it approaches convergence and as the thermody- 
namic limit is taken. Central to this approach is the manifestation of extensivity through a 
description based on the Cumulant Generating Function, which we take to be given. We then 
derive a set of general integral equations which define the scaled Lanczos functions in the 
thermodynamic limit, which can be explicitly and exactly solved for certain integrable mod- 
els, or employed in a truncated manner for non-integrable models. An alternative formulation 
is also given which expresses the equivalence of the Lanczos Process with the continuum Toda 
Lattice Model treated as a boundary value problem. Finally we state some general results 
concerning the behaviour of the Lanczos functions. 
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II. The Lanczos Process, Orthogonal Polynomials and Mo- 
ments 



The Lanczos Algorithm or Process]^, ||[ U begins with a trial state l^o) appropriate to the 
model and the symmetries of the phase being investigated. From this the Lanczos recurrence 
generates a sequence of orthonormal states {|V>n}}$£Li and Lanczos coefficients {a n (N)}^ =0 



and {P n (N)}™ =l , thus 



H\lp n ) = (3n\lpn-l) + a n \lp n ) + f3 n+ i\lp n+1 ) 

with the Lanczos coefficients being defined 

a n = (ip n \ft\il>n) , 

Pn = (lp n -i\H\l/j n ) . 



(1) 



(2) 



We distinguish a total or extensive operator or variable such as H from its density or intensive 
counterpart by H . In this basis the transformed Hamiltonian takes the following tridiagonal 
form 



( a fi\ 

Pi Oil P 2 



\ 



\ 



Pn-1 "n-l Pn 

Pn OL n ) 



(3) 



As such the Lanczos process is one of the Krylov subspace methods 0, in that at a finite step 
n, the eigenvectors belong to the Krylov Subspace Span{|^o)> H\ipo), H 2 \ipo), . . . H n \vpo)}. 

In the Many-Body context one would iterate the Lanczos Process until termination 
whereupon the Hilbert space is exhausted (at this point one of the P nT = 0, where tit is the 
dimension of the Hilbert space in the sector defined by the ground state), or until the process 
has converged according to some arbitrary criteria n — > tiq- Then one would perform the 
thermodynamic limit — > oo where it should be understood that the above conclusion of 
the Lanczos process is also dependent on the system size, that is to say riT{N),nc{N). These 
cutoffs are monotonically increasing functions of the system size so they will all tend to oo 
in the thermodynamic limit as well. Taking the limits in the reverse order clearly leads to 
nonsensical results, as taking N — » oo with n fixed produces a n — > c\ and p n — > 0. The great 
virtue of the Lanczos process is that it can be shown to converge essentially exponentially 
fast with respect to iteration number, using the Kaniel-Paige-Saad exact bounds |], 0, ||] for 
the rate of convergence. This means that convergence occurs within a very small subspace 
of the total Hilbert space, so that uq « tit- 

The Lanczos process is entirely equivalent to the 3-term recurrence for an Orthogonal 
Polynomial System [l0|, 11], however we consider a slight generalisation of the preceding 



process to one with a single parameter evolution (a "time" t). In this construction we are 



continuing a development begun by Lindsay [12] and Chen and Ismail] 13], which will lead to 
some powerful tools in treating the Lanczos process. The measure, or that component which 
is absolutely continuous, is defined by the weight function 



w(e,t) 



-u(e)+tNe 



(4) 



on the real line e£l. Our system under study is described by the initial value of the system at 
t = and often we will suppress this argument for the sake of simplicity. This measure defines 
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a system of monic Orthogonal Polynomials {P n (e,t)}^ =0 with an orthogonality relation 

/+oo 
de w(e,t)P m (e,t)P n (e,t) = h n (t)S mn , (5) 
-oo 

and normalisation h n (t). This is equivalent to the following three-term recurrence relation 

P n+1 (e,t) = (e - a n (t))P n (e,t) - p 2 n {t)P n ^(e,t) , (6) 

with the recursion coefficients a n (t) real for n > and 0^(t) real and positive for n > 0. By 
convention we take @l = 1. It can be readily shown that the Lanczos coefficients are given 
in terms of the normalisation thus 

an{t) = wmjt K{t) ' m 

K{t) { ' 



K-\{t) 



The direct connection between the Lanczos Process and the OPS are given by the determinant 
relation of the characteristic polynomial 

P n+1 (e) = (-)" +1 |r n -e/ n+ i| , (8) 

so that the zeros of the Orthogonal Polynomial are eigenvalues of Hamiltonian. 

Some comments are in order regarding the differences, or more accurately the special 
character, of these Orthogonal Polynomials with respect to the generic OPS or with some of 
the scaling versions of OPS|Q. These OPS have been termed Many-Body OPS, but could 
be equally described as extensive OPS. They all have an additional, essential parameter to 
the generic OPS, the system size N, which appears in both the gross scaling factors (the 
'external' scaling such as in the energy densities e defined by E = Ne), but also internally 
in the 3-term recurrence coefficients, in the Polynomials themselves and in other derived 
quantities. The internal dependence in the Lanczos coefficients on the system size is not at 
all apparent and the most transparent way that extensive scaling properties can be exhibited 
is through the Cumulant Generating Function, (CGF), which hitherto has played no role 
in Orthogonal Polynomial Theory. In fact the CGF is central to this class of OPS rather 
than the moments, and is in a practical sense the starting point in any application of the 
Formalism to physical Models. For all models it is clear that the ground state energy Eq 
is proportional to N and unbounded in the thermodynamic limit, and similarly the total 
Lanczos coefficients (as opposed to the densities) are unbounded as n — > oo for fixed N. 
When everything is recast in terms of densities the spectrum is bounded below by eo and in 
many models will also be bounded above, and similarly the density Lanczos coefficients are 
bounded. Another difference that Many-Body OPS exhibit in comparison to general OPS is, 
as we have noted above, the three-term recurrence will terminate exactly at n = tit, although 
this will never present any problems as this is exponentially large. 



The Lanczos process is intimately connected with the Hamburger moment problem 1 15, 
16 1, via the Resolvent operator 



R(e) = (-^-r) e i Supp[dp] . (9) 
e — H 

Its formal Laurent series establishes a direct link with Hamiltonian moments 

oo 

Mi 



^) = £-St> (io) 



e 
i=o 
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where these moments are defined as expectation values with respect to the trial state referred 
to above 



Mr 



(H n ) 



Mo 



1 



The resolvent has a real Jacob i-fr action continued fraction representation! 17, 18 1 

!3 2 



R(e) 



J:s -n=0 



a r 



(11) 



(12) 



with elements coming from the Lanczos coefficients. 

An equivalent description to that of the Hamiltonian moments is to formulate everything 
in terms of cumulants or connected moments^, [H]] {^n}^i> an d to ignore all corrections 
which vanish in the thermodynamic limit N — > oo. Cumulants scale directly with the size of 
the system so that for the extensive Many-Body Problem we have 



in the ground state sector, or 



v n = c n N + o(l) 



c n N + m n + o(l) 



(13) 



(14) 



in any other sector [21]. This also means that no finite-size scaling can be performed given 
that only the limiting quantities are retained here and boundary condition effects do not 
appear. The foundation ingredient is the Moment Generating Function which is related to 
the Cumulant Generating Function in the following way. 

Definition 1 The Moment Generating Function (MGF) M(t) and the Cumulant Generating 
Functions ( CGF) F(t) are defined by, 



M{t) = (e 



tH\ 



n=0 



IV. 



exp 



\n=l 



exp(NF(t)) . 



(15) 



Some examples of Cumulant Generating Functions include the isotropic XY model using the 
z-polarised Neel state as the trial state [22] 

1 /-w/2 



K Jo 



dq log cosh(t cos i 



(16) 



and the Ising model in a transverse field using the disordered state as the trial state, and 
coupling constant ar[23| 



W 2jt Jo 



dq In 



cosh(2te^ 



(cos q+x) 



sinh(2te 9 



(17) 



where the quasiparticle energies are defined by t 2 a = 1 + x 2 + 2x cos q. 



Definition 2 The Determinants of the Moment Matrices A n (£) for n > are defined by the 
Hankel form - 



A n (t) = \M^- 2 )(t)\^ 1 . 



(18) 
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The direct relationship from moments to the Lanczos coefficients which is established in this 
way is via the construction of a sequence of Hankel determinants of the Moment Matrices 
and their Selberg-type integral representation || 



j r+oc n +l 

V /■ J-oo k=1 Ki<j<n+1 



%<]<_ 

These determinants are related to the normalisations via 



Definition 3 Our final definition, that of the Lanczos L-function, is 



^ = ^gl, (2!) 



for n > 1 and L (t) = M(t). 
The converse result is then 



k=0 

for n > 1. From these the Lanczos coefficients are given simply by 

3 



(3 2 n {t) = L n (t) 



(19) 



A n (i) = Y[hj(t) . (20) 

j<n 



A n (t) = N n ^ J] L n k +1 ~ k (t) , (22) 



(23) 



Theorem 1 T/ie equation of motion for the Lanczos L-functions is 

M*) = ^E^A 2 iogW*) • ( 24 ) 

with the initial condition on the recurrence given by logLo(i) = NF(t) for all t. 

The advantage of introducing evolution into the Lanczos Process is that Sylvester's Theorem 



applied to the Hankel determinants [24] 



A n+1 (t)A n _ x (t) = A n (t)K(t) - {A' n (t)) 2 (25) 



so that the theorem follows directly from this.D 

The first few members of the Lanczos -L-sequence are 



L 1 (t) = iF"(t) , 



JV W ' AT2 (i?(2)' 



2 



The consequence of Sylvester's theorem for the evolution of the A n is the following 
theorem 
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Theorem 2 The A n (t) obey the following differential- difference equation 

exp {log A n+1 + log A n _x - 2 log A n } = D 2 t log A n , (27) 

with the boundary value logAo = NF(t) and conventionally A_i = 1. 

This follows directly from Sylvester's Identity. □ 

This evolution equation is just the finite Toda Lattice equation of motion |25|| , and this 
point has been previously noted in Ref. [13|. 
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III. Scaling in the Thermodynamic Limit 



As was discussed earlier there are two limiting processes that one must consider when the 
thermodynamic limit is taken in the Lanczos Algorithm, both n, N — > oo, and the issue then 
is what mutual relationship exists between them in the limit. One can view this limiting 
process in the 1/n vs. 1/N plane and then consider along what types of paths must one 
approach the origin. We shall find that the general relationship is n, N — > oo with s = n/N 
fixed, although for systems at criticality it seems inevitable that s will become unbounded 
in the analysis. A consequence of these ideas is the confluence property of the Lanczos 
coefficients as n, N — » oo at fixed s = n/N 



a n (N) = a(s) + O(l/A0 , 
(3 2 n (N) = f3 2 (s) + 0(l/N) . 



(28) 



There are a number of ways to see this approach to the thermodynamic limit. 

Using the explicit forms connecting cumulants and moments, and a direct evaluation of 
the Hankel determinants one can prove [2^] for general n and iV that the Lanczos coefficients 
have a leading order scaling in s = n/N for the first two orders of an expansion in large N. 
Actually this expansion is valid for all n not just for large values and thus includes all the 
subdominant contributions. Thus 



a, 



--ciN + n 



+ -n(n-l) 
+ ... , 



3c|-4c 2 C 3 C4 + C2C 5 

24 



f 

N 



(29) 



for n > 0, and 



0l 



nc2N + —n(n 



H — n(n 
6 

+ ... 



l)(n 



c 2 c 4 - 



,,2l 



- 1 2c| + 2 1 C2 C3 C4 — C4 — 6C2 C3 C5 + C% C6 

2^ 



f 

N 



(30) 



for n > 1. However this approach cannot be generalised to higher orders and therefore for 
the full exact Lanczos coefficients. The first two terms in the above expansions were also 
proven by Lindsay using the Sylvester Identity in the statistical context |i~2]] but no further, 



while this form for the higher terms (but finite numbers) was conjectured in Reference [27]. 
We shall find that use of the Sylvester Identity allows one to very easily recover this result, 
to in fact go to much higher orders in constructing explicit forms and to prove this type of 
scaling in a completely general way. 



Lemma 1 The Lanczos L-function L n (t,N) is a rational function of 1/N for fixed n, and 
all t. 



The Difference-Differential Eq. 
Also for fixed n we have 



24| ) is of finite order in j/N and t, so the result follows. □ 



N 



lim L n (t,N) = 



(31) 
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and specifically the leading order term is 0(N 1 ) which arises from the j = n term in the 
sum. Therefore we can expand this function in a descending series in N' 1 , thus 

M*,A0 = £%^> (32) 
p>l 

and defining the connected series related by 

P >1 \ P >i / 

This last relation can be rendered into an explicit form 



«•'»■ - E E«-D'IIp( z tf i )"- < 34 > 



It is actually necessary to perform an expansion of this type because it combines the iteration 
number (n) dependence of the numerator and denominator which are both essential in the 
following results. 

Then one can establish a hierarchy of equations for these coefficients 
l nl (t) = nF"(t) , 

n-l 

ln2{t) = Y j jD 2 t \ogl n . jl {t) , 

3=1 (35) 
n-l 

l n P {t) = jm'^_ jp _ 2 (t) for p > 3 , 

3=1 

for n > 1 whilst for n = we have / n p(i) = as Lo(t, N) = exp(NF(t)). The first members 
of this hierarchy can be easily solved for yielding 

l nl (t)=nF"(t) , 

1 i?(2) F (4) _ rp(3))2 
ln2(t) = -n{n- 1) , 

2 FT (36) 

^3(t) = T2 n{n - l)(n - 2) f (i?(2)) l 3 j , 

and from these it is easy to establish the leading order terms already found in Eq. (2S,30p. 

Lemma 2 The hierarchy coefficients l np (t) , m np (t) are polynomials in n. 

These coefficients are constructed from a finite difference equation in n. □ 

Theorem 3 The hierarchy coefficients l np (t) , m np (t) are polynomials of degree p in n. 

This is proved by induction on p using the hierarchy equations. If we take lj q (t) to be of 
degree q < p — 2 in n then similarly for mj q (t) and m" (t). Now for any polynomial P(n) of 
degree p — 2 in its argument then 

n-l 

E-? p ( n --?) > ( 37 ) 

3=1 
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is apth degree polynomial. Thus the recurrence, Eq. (j35|), establishes that l n +ip is also apth 
degree polynomial. □ 

From this result it is clear that the limiting forms of the Lanczos coefficients a n (N),f3^(N) 
exist when n, N — > oo with n/N fixed. If the ratio is not kept constant in this limiting op- 
eration, say with n = o(N) then the Lanczos coefficients will vanish in the limit, while if the 
reverse is true N = o(n) then there will be divergent terms in the limit. 

Given that the scaling Lanczos coefficients have been established then all the exact 



theorems for the ground state properties^, 29] that were predicated on this result now are 



established. The first example of these theorems was the one for the ground state Energy 
Density, 

e = inf ■ \a(s) - 2(3(s)\ , (38) 

which also has an analogue for the top of the spectrum, if this exists 

eoo = sup [a(s) + 20(a)] . (39) 

For many models these Lanczos Functions will be bounded on the positive real axis, and have 
limits as s — > oo on the real line. So there is a superficial similarity to classes of Orthogonal 
Polynomials whose 3-term recurrence coefficients have limiting values, such as the S-class, 
the M-class, or the M(a,b) classes|14]. 
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IV. The extensive Measure 



It is necessary to determine the OPS measure, its weight function w(e), and this is not 
generally known at the outset, but rather the Cumulant Generating Function is. In fact it 
seems to be the case that the measures are not exactly expressible in simple terms, but the 
CGF or characteristic functions are. There is of course a direct route from a model system 
and a trial state to the Lanczos coefficients, but from many points of view including practical 
considerations the route beginning with a cumulant description is more useful. 

Theorem 4 Given that the cumulant generating function F(—t) is analytic for 3t(t) > 
and in the neighbourhood of the origin t = the OPS weight function w(e) has the following 
asymptotic development in the thermodynamic limit N — ► oo ; 



tt ' (e) = Vi^ e " H!+F,£,l+0(iV " /2) ' (40) 

where the function £(e) is defined implicitly by 

e = F'(0 . (41) 
Starting with the definition of the cumulant generating function F(t) 

{e tH } = eM NF(t)} = exp jjvf; ^t" j . (42) 

We assume here that this infinite series is not just formal but actually exists, that is it has 
a finite radius of convergence in addition to its analytic character for 5?(t) < 0. However 
the Moment Generating Function is simply the analytic continuation of the characteristic 
function and this continuation is possible given its analyticity, so that a Fourier inversion of 
this will yield the weight function, 

TV ril+oc 



AT rt^+co 

w(e) = / dt e N [-W(.t)] , 

Jil ~Z (43) 

AT />7+JOO 

= A / dt e mte+F(-t)] K(7) > 



One does not require the exact inversion but only the leading order in N in a steepest descent 
approximation. In an asymptotic analysis the relevant function is 

g(t) =te + F(-t) , (44) 

which is analytic for all 3t(t) > 0. We will assume the existence of a stationary point which 
occurs at to 

e = F'(-to) , (45) 

and is assumed to be unique. This point is evidently real because the energy density is real 
and the CGF is a real function of a real argument (here we define £ = —to for convenience). 
One requires the inversion of this relation for £(e) and this is guaranteed by the Implicit 
Function Theorem because F^ 2 ^(£) > 0. This latter condition also implies that the saddle 
point is of order unity. Indeed one clearly has the case of F^ (t) > for real values of t in 
the neighbourhood of the saddle point and F^ (t) < for imaginary values of t in the same 
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neighbourhood. Thus the path of steepest descent through the saddle point is parallel to the 
imaginary axis. One can then apply the standard saddle point analysis, see Wong[j30|l Section 
II. 4, to arrive at the stated result. □ 

The corresponding example of the saddle point equation for the isotropic XY model is 



1 r*l* 

e = — dq cos q tanh(£ cos q) , (46) 

JO 

and that for the Ising model in a transverse field is 



1 r ^±g^+tanh(2^ 9 ) 

6 = - / d ^ h^Ta ' ( 47 ) 

it J 1+ + eg tanh(2£e g ) 

The first of the more obvious properties concerns the convexity of the measure arising 
in the thermodynamic limit, 

Theorem 5 The leading order of the negative logarithm of the weight function u(e) is convex 
for all real energies e. 

This follows from the relationship of u(e) to the stationary point 

ju(e) = Nfa) , (48) 

and the definition 

e = F>(0 . (49) 

Now it can be easily seen that F"{t) > for t real and the Hermitian Hamiltonian using the 
definition of F(t) in terms of the expectation value NF(t) = ln(exp(ti?)).D 

Some detailed, yet general information, concerning the extensive measure in the neigh- 
bourhood of the ground state is available. This arises from consideration of the overlap of 
the trial state with the true ground state[31], and its relation to the Horn-Weinstein function 
E{t) = F'(-t) via 

\(*GSm\ 2 = ™vl-N J™ dt [E(t) - E(po)]} . (50) 

In general the limit E(t) as t — > oo will exist, and is the ground state energy, and so the 
asymptotic properties of E(t) for $t(t) > as this tends to infinity is a means of classifying 
systems. This equivalent to the asymptotic properties of e(£) — e(— oo) as £ — > — oo (we 
denote the Ground State Energy by eo, which is also the same as e(— oo)). In general the 
overlap is non-zero, so that E(t) — E(oo) G ^ 1 [0, oo) but it is possible at isolated points that 
this is not true (critical points in the model for example) and the overlap may vanish. For 
example the overlap squared in the case of the isotropic XY model is 2~ N / 2 and that for the 
Ising model in a transverse field is 



2tt J \ 2e q 



Where the overlap is non-zero then several possibilities for the asymptotic behaviour exist, 
which do actually arise in the exact solutions of the example models - 
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• gapless case, isotropic XY and critical Ising Model in a transverse Field, Ref.|3^, 23]: 
At a critical point, the first excited state gap vanishes and 

e-eo~yl|£r 7 , (52) 

as £ —* — oo and if the overlap is finite then 9?(7) > 1. Therefore the weight function 
at the bottom of the spectrum takes the following form 

w(e) ~ (e - eo)"^ e xp |iV-A-( e - co) 1 " 1 /? J , (53) 

This measure is integrable on (eo, eoo) because of the above condition $£(7) > 1 and has 
a branch point at the ground state energy eo- 

• gapped case 1, Ising Model in a transverse Field, in the ordered phase with the disor- 
dered trial state, Ref. p3fl : 

if the gap is finite then one possibility is that 

e - e ~ Ae~ m , (54) 

as £ — > —00 and where the excited state gap A > 0. One can show that the weight 
function near the bottom edge of the spectrum is analytic having the form 

w(e) =_ . (55) 

r(jV^l + l) 

• gapped case 2, Ising Model in a transverse Field, in the disordered phase with the 
disordered trial state, Ref.[23|: 

and yet another type of gap behaviour exists 

e-e ~ A\£\-~<e- m , (56) 
The leading order behaviour of the weight function in this case is 

W(e) ~ (6 - eo )-V2-7V(e-e„) { _ log{e _ ^-N^e-eo) > (57) 

which again has a branch point at the bottom edge of the spectrum. 

So generally we find the support of the measure is bounded which excludes a number 
of weight function types such as the Freud or Erdos weights, but that the weight functions 
belong to the Szego class on [eo,eoo]j 

3 , log We) 
de s v I = > -00 . (58) 

V[eoo - e][e - e ] 
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V. Exactly Solvable Lanczos Process 



In this section we derive how the exact Lanczos functions a(s) and f3 2 (s) can be constructed 
directly from the knowledge of the connected Moments or Cumulants, or more specifically 
from the Cumulant Generating Function. This is the initial data that one uses in any 
analysis of quantum Many-Body Systems with this approach, and for soluble models the full 
Generating Function may be available. However if this is not the case then one would use a 
set of low order Cumulants, up to a given order. 

As a first step we recast the Hankel determinants into Selberg Integral form, from the 
classical result [01 



A„(t) 



1 



(n+1)! 



n+1 



(59) 



' k=l l<i<j<n+l 

For the steps leading to the two conditions which will define the Lanczos functions we 
follow Chen and Ismail[13[. A similar approach, but just confined to the evaluation the 
Hankel determinants, was taken in References [33, 34]. The Hankel determinant can be 
recast into the form of a partition function, which is, 



A n (i) 



1 



(n+1)! 



n+1 



n+1 



n+1 



n+1 



exp 



Y u(ei) + Nt Y ei + 2 Y ln l e * - e i 

i i i<j 



(60) 



One should observe that both Yl'i u ( e i) an d NtY^i € i are of order (n + 1) AT whilst the 
remaining term in the argument Y^i<j m l e « — e j \ * s °f or der (n+1) 2 , so that the only relative 
scaling that remains nontrivial is one in which n/N is fixed. The alternatives would lead 
to completely trivial consequences. The leading order term for this Hankel determinant as 
n, N — > oo is given by a steepest descent approximation (see Ref. [3C] section IX. 5) 

-1/2 



A n (t) 



(27T) 



n+1 



(n+1)! 

where the function /(e) is defined as 



de°de 



o 



1 + 0(l/n, 1/JV)] , 



n+1 



n+1 



n+1 



NtY^i 



i i 

and the saddle points {e?}^ are given by 

u'{£) =Nt + 2 



\7I ) 



(61) 



(62) 



i<j 



n+1 



(63) 



One can easily show that the Hessian in Eq. ( |6l"D is positive definite given that u(e) is convex. 
One can carry the continuum limit further by describing the saddle points as a charged fluid 
whose dynamics are governed by an Energy Functional F[a] 



exp 



n,N—>oo 



exp(-F[cio]) , 



(64) 



with a charge density <r(e) defined on an interval of integration which is to be determined, 
/ = (e_, e + ). The energy functional takes the following form 



F[a] = J dea(e) [u(e) - Nte] - J^de J de'a(e) ln [e - e'\a{e') , 



(65) 
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where the single particle confining potential is controlled by the OPS measure and the two- 
body interaction is a logarithmic type. The result of minimising this Functional yields the 
following singular integral equation for the Charge Density 

u'(e)-Nt = 2PV fde'^^ . (66) 
J I e — e 

The solution of this integral equation for the Minimal Charge Density 00(e) can be found 
exactly and is 

= sMEWEEi PV \x w- m . (67 ) 

There are two conditions arising from this solution - 

• the first is a Supplementary Condition which is necessary for the charge density solution 
to be well defined throughout the interval / 

J I V( e + - £ )( e - e -) 

• and the Normalisation Condition which simply counts the number of Lanczos steps 

1 f , u'(e) -Nt 

n = — de e — w . 69) 

Using this solution for the charge density one can substitute this into the original defining 
equations for the Hankel determinants (the leading order approximations) and establish that 
the Lanczos functions are simply defined by the interval I in this way, e± = a ± 2/3. 

Theorem 6 The Lanczos functions are given implicitly by the two integral equations 

= Ja-20 de ^/4f3 2 -(e-^ ' (70) 

s = — de ; gl ' = , (71) 

where the model dependent equation for the stationary point £(e) is given by Eq. ^4W- 



This theorem follows from the previous conditions, namely Eqs. (pq, 39 ) , and the result for 
the logarithmic derivative of the weight function, 

u'(e) = iV£(e) + 0(logiV) . (72) 

□ 

Usually this later equation for the saddle point is also an implicit equation and invariably 
a nonlinear one. In our derivation the scaling s = n/N remains finite whilst n, N — > oo 
emerges naturally and in fact it is difficult to see how one could avoid this confluence. 

We now give an alternative result for the Lanczos functions which is based on the time 
evolution of the Lanczos L-function. 
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Theorem 7 The Lanczos L- function, in the thermodynamic limit is the solution of the fol- 
lowing integro- differential equation 

L(s,t)= I dr rD 2 log L(s -r,t) + sF {2 \t) , (73) 
Jo 

and the two Lanczos functions are derivable from this via 

a(s) = [ dr D t logL(r,0)+F'(0) , 

Jo (74) 

P 2 (s) = L(s,0) . 

The integro-differential equation is simply derived from the discrete recurrence, namely 
Eq. (|24]), after making the observation that the j = n term involving Lo(t) has to be sepa- 
rated from the sum because it encompasses the initial conditions and is itself not generated 
by the recurrence. □ 

Finally we give a result equivalent to the theorem above, but which involves only scaled 
forms of the Hankel determinants A n (N,t) and is the differential analogue of the above 
Theorem. 

Definition 4 We make the following definition for S(n, N, t) in terms of the Hankel Deter- 
minant, 

A n (N, t) = N n ^ [5(n, N, t)] N * , (75) 
for n > 1 and A (t) = [5(0,*)]^. 

Lemma 3 The function 5(n,N,t) is well defined in the scaling limit n, N — > oo. 

This follows naturally from the relation of the A n (t) and the Lanczos L- function as given in 



Eq. ( 22 ) , and the well defined scaling of this latter function as demonstrated in the Theorem 
3 above. □ 

Then we have the following result - 

Theorem 8 The Lanczos 5 (s,t)- function satisfies the following partial differential equation 
in the thermodynamic limit 

exp{L» s 2 log S(s,t)} = Dj log 5(s,t) , (76) 

with the boundary condition 

lim iog5M) = VfGM+) 

s^0+ S 

The Lanczos functions are given by 

a(s) = D t D s logS(s,t)\ t=0 , 
P 2 (s) =exp{D s 2 log5(s,0)} . 



Using the scaling relation above, Eq. (|75|), and the equation of motion for A n (i), Eq. 
the result follows. □ 

These last two theorems relate to the dynamics of a nonlinear continuum Toda Lattice 
in one space domain s E M + and one time domain t, with boundary conditions defined at the 
origin s = for all times t by the cumulant generating function F(t). The object is then to 
find the Lanczos functions a(s), (3 2 (s) from a solution of this system, wherein these functions 
are directly related to the solution at a given time t = over all spatial points s > 0. 
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VI. The Taylor Series Expansion 



The investigation of the Taylor series expansion of the Lanczos coefficients about s = 0, is an 
essential element in the application of this Lanczos method, as was indicated earlier, where 
one has only a finite set of low order cumulants available, say for non-integrable models. 
Therefore in this case one can only construct a truncated Taylor series expansion and so 
issues concerning convergence, the radius of convergence of the series, and whether one can 
extrapolate immediately arise. In addition one would like a direct algorithm relating the 
cumulants to the Lanczos functions from a purely practical point of view. 

We define the Taylor series expansion of the two Lanczos functions by two new sequences 
of coefficients, 



a(s) = ci + J2 a nS n+1 , 

(79) 



n=0 



In order to find these coefficients one could use either of the two general solutions for the 
Lanczos process, Eqs. (70, 7^) or Eq. (|73|), and the two methods are presented below. 



The first step involves the inversion of the following Taylor series expansion 

e = ci+ j2^fe, (so) 



n=l 

for £(e), namely the coefficients e k appearing in 

t = ^2e k (e- Cl ) k . (81) 



k=l 



The coefficients c n appearing in Eq. 80 are the cumulant coefficients. The existence of this 
inverse function is guaranteed because the second cumulant C2 > in all systems and we 
assume that the saddle point function, Eq. (|49"|), is analytic in the neighbourhood of £ = 0. 



The next step involves the solution of the two recurrences 

[k/2] 



E e * E (? )^(«-ci) fc - 2m (4/3 2 r , 

^— ' \2mJ ml 

/ i :n I I •• (82) 



L(fc-1)/2J f , . 

2s = Ye k Y ( k ) \^h±l (a - c 1 ) fc - 2m - 1 (4/3 2 r+ 1 , 
^ n \2m+l (m+1 ! V ' V H ' 

k=l m=0 x / v ' 

which are used to solve for the coefficients a n , b n appearing in Eq. (|79|) . 

In the second method we define a continuum version of the coefficients that are defined 
in Eq. (^) in the following way 

log inx = lo q 1 + E lj fA - E m ^) sP . (83) 

and the inverse of Eq. (^) in an explicit form 

¥= E Up<- < 84 > 

l l ^ ,• HI- 



17 



From these relations one can find a hierarchy of equations for these coefficients 



h(t) = F"(t) , 

i?(2) i 7(4) _ (_p(3))2 



hit) 



2(F(2))2 ' (85) 



mi 1 



Thus one can verify from the solution for the initial value problem above that the general 
Taylor series coefficients are given by 

2n+l 

Ah2n+1 i=0 , 

2n ^ ' 



[(n+l)!] 2 c^a n = Y, A ^ X ) II c % 

Ah2n+1 i=0 
2n 

{n+l)\n\cf- l b n = B(n;X)Hc a 2 ^ , 



Ah2n i=0 

where the coefficients labeled by the partition A = (l ai .2 a2 . . . i ai ), denoted by A(n; \),B(n; A), 
are listed in Table(||) of the Appendix. There are constraints operating in the above equations, 
namely X^i" 1 za i = X^=o 1 a * = 2^+1 for the first relation and Y^i=i = Si=o a * = ^ n ^ or 
the second. 

Clearly the Taylor series expansion of the Lanczos functions has low order coefficients 
which are constructed from the low order cumulants, and is a form of a linked cluster expan- 
sion. However it is not just a simple linked cluster expansion as in the Taylor series expansion 
of the Cumulant Generating Function, but involves a subtle interplay and cancellation of all 
cumulants below a given order. 
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VII. General Properties 



There are some very general properties that the Lanczos process in the thermodynamic limit 
and the associated Lanczos functions satisfy and we examine these now. Some are quite 
obvious and not particularly surprising, however we state these for completeness sake, while 
there are some other properties which are not so immediate but very important nevertheless. 

The next, and natural, property concerns the monotonicity of the two envelope functions 
e±(a) =a(s)±2(3(s). 

Theorem 9 The envelope functions e + (s), e_(s) are monotonically increasing and decreasing 
functions of real, positive s respectively. 



This follows from a recasting of the normalisation condition in the following way 



2ns = j <+ deV[e(£+)-e(0][e(0-e(^)] , (87) 



where the £± are defined by e(£±) = e±. Now it is straight forward to write the explicit forms 
for the derivatives of the envelope functions with respect to s as 




<t)-<t-) ' 



so that the stated properties are evident. □ 

It is clear that the envelope functions e±(s) are bounded in the following ways, e_(s) > eo 
and e+(s) < e^. 

The 3-term recurrence which serves as one of the definitions of the Orthogonal Polyno- 
mials themselves is now going to take a definite limiting form when n, N — > oo such that s is 
finite. This is going to lead to a scaling form for one set of the Polynomials themselves, which 
would be more correctly termed orthogonal functions p(s,e). Heuristically one can see how 
this arises by the following argument. If one ensures that Lanczos densities are employed and 
the following scaling of the polynomials thus P n (E) = N n p n (e), then the 3-term recurrence 
becomes 

Pn+i(e)/p n (e) + 0i — -— p -— = e-a n . (89) 

Pn(eJ/p n _i(eJ 

Now these ratios are approximated by 

M(i| -=p(^Mv)), (90) 



for arguments e € C\Supp[dp]. So that in the asymptotic regime the recurrence becomes 
19, , A 2 fid 



CX P( ^^ ln ^ s ' e ) ) + ^ ( s ) ex P ( -^7^ ln ^( s ' e ) ) ~e-a(a) , (91) 



whose solutions are 



p ± (s,e) ~p(0)exp|AT J dtln^ e - a(t) ± y/(e - a(t)) 2 - A(5 2 {t) | . (92) 
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These are the corresponding results for the ratio P n (x) / P n +i(x) or n-th root y/P n (x) asymp- 
totics of generic Orthogonal Polynomials as n — » oo|| [ll], |3^, ^6, 37], or the scaled Orthogonal 
Polynomials [S3] , but are rather different due to the particular nature of Many-Body Orthog- 
onal Polynomials. 

Theorem 10 Given the scaling behaviour of the Lanczos coefficients, and that they are 
bounded for n,N — > oo, then the n-th root of the denominator Orthogonal Polynomials p n (e) 
have the limiting form uniformly for e in compact subsets of C\Supp[dp] . 



p(s,e)= lim \p n (N,e)\V N 

n,N—>oo 



exp 



l r 
dt In- 
2 



-a(t) + 



(93) 



The proof of this parallels the one constructed by van Assche in Ref. 
Turan Determinants, 



One can show that these obey the following recurrence relation 

D„ = $^ D n-l + {a n -a n -i)p n p n -i + (0%- Pn-t)p n Pn- 



through the use of 
(94) 

2 (95) 



Using the partial fraction decomposition of the ratio of two successive Orthogonal Polyno- 
mials one can also find a bound on this ratio 



Pn-l(<0 



Pn(e) 



< — V n 
d 



(96) 



for all e G if where the compact set if C C\Supp[dp] and d is the distance between this set 
and the interval [eo, £oo]> and C is a positive constant. Using Eq. (|95|) we have 



c 2 



D 



n-l 



Pl-1 



(97) 



Given the scaling form of the Lanczos coefficients the ratio \D n /p^\ —* as n, iV — > oo 
uniformly in e whenever d is large enough. This means that \p n _i/p n \ and |pra/Pn+i I tend to 
the same accumulation point which we denote by p(s, e). This point is given by the solution 
of the quadratic equation p + j3 2 {s)/p = e — a(s), and the positive branch of the solution 
must be taken as p — > oo when e — > oo. The functions p(s, e) are analytic functions of e £ K 
which are uniformly bounded, so the restriction on d can be lifted to being only non-zero. 
The behaviour of the n-th ratio then gives the n-th root behaviour directly as 



\p> 



il/JV 



exp 



k=l 



Pk{£) 



Pfc-i(e) 



(98) 



The asymptotic behaviour that we have found applies to the denominator OP only as can be 
seen from the observation that p\ = e — c\ and p2 = (e — c\) 2 — c^/c2N(e — c%) — 02/N, while 



r-a(s) + V(e-a(s)) 2 -4P 2 (s) 
This establishes the result. □ 



1 



+0 e - c\ 



{(e-a) 2 -c 3 /c 2 N(e-c 1 )-c 2 /N) . (99) 
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VIII. Summary 



In this work we have demonstrated the general scaling behaviour of the Lanczos Process as 
applied to Many-Body Systems when the process is taken to convergence and the thermody- 
namic limit taken. We also find explicit constructions of the limiting Lanczos coefficients in 
two equivalent formulations, from an initial exact solution of the moment problem, that is to 
say the cumulant generating function for the system. There are explicit examples where the 
CGF can be found and the whole Lanczos process explicitly realised. Furthermore we have 
given the corresponding results for the associated Orthogonal Polynomial system and the 
measure in this regime, quite generally. However we must emphasise that these results apply 
only to the bulk properties, that is to say the ground state properties that scale extensively 
and the spectral properties in the interior (the "bulk") of the spectrum. So this does not 
include the delicate scaling behaviour at the edges of the spectrum, nor in the neighbour- 
hood of singularities - this theory would have to be extended to treat the excited state gaps 
near the bottom of the spectrum. A number of general theorems are given which constrain 
the behaviour of the Lanczos functions, and the process in general. We also indicate how 
a number of such constraints operating can lead to some concrete realisations or scenarios 
that the Lanczos process can present, namely its behaviour at a critical point in the model 
under study. This is a significant step on the way to the goal of a rigorous classification of 
Many-Body Systems in terms of their character via the Lanczos process. Other important 
questions that arise in the treatment of non-integrable models, for which the general results 
presented here have suggested some answers, are the questions of the choice of trial state, 
the rate of convergence of the truncated Lanczos process and how one might accelerate its 
convergence given some independent qualitative knowledge. 
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Appendix 

We list here the coefficients of the Taylor series expansion for the Lanczos Coefficients, labelled 
by the partitions of integers, according to the definition of Eq.(86). 
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Table 1: The coefficients in the Taylor series expansion for the Lanczos functions a(s) and 
/9 2 (s), as defined in Eq.(p6j), and the labels denoting the partitions A of the positive integers. 
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